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Abstract 

In this paper, a second order wage equation is developed and solved by the method of variation 
of parameters. The subsequent wage function is then analyzed and interpreted for stability. 
Speculative parameters, which operate freely dictating employers' expectations, are included in 
modeling this equation. The variation of these parameters causes both stability and instability of 
the wage function depending on circumstances. Where the wage function is exponential, 
asymptotic stability towards the equilibrium wage rate is observed but where it consists of both 
exponential and periodic factors, the time path shows periodic fluctuations with successive 
cycles giving smaller amplitudes until the ripples die naturally. It has been realized that where 
the wage rate is determined by free market forces of demand and supply, volatility in wage rate 
may be observed if not controlled. This may increase uncertainties and cause anxiety about 
investment and employment in the economy. The paper therefore proposes government 
intervention by creating a middle path in which wage rate is allowed to oscillate freely within a 
narrow band managed by employers in consultation with the workers under the watch of the 
government. 
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1. Introduction 

Wages mean the reward for labor services. In [6] wage is defined as a fixed regular payment 
earned for work or services, typically paid on a daily or weekly basis. In [9] wage is viewed as 
payment for labour services to a worker, especially remuneration on an hourly, daily, weekly or 
by piece. It also views wage as a portion of national product that represents aggregate paid for all 
contributing labour and services as distinguished from the portion retained by management or 
reinvested in capital goods. In [10] wage is defined according to wages act of 1986. It is the sum 
payable to an employee by an employer in connection with that employment. It includes fees, 
bonuses, commissions, holiday pay or other emolument relevant to the employment whether 
specified in the contract of employment or not. 

In this paper, we consider modeling a second order differential wage equation. For example, 
considering deterministic price adjustment model in [7] fixed supply and demand functions at 
instantaneous price for security is discussed. It is argued that at equilibrium asset price, quantity 
demanded equals quantity supplied. This is discussed using fixed demand and supply curves 
while price is kept constant. It also asserts that away from the equilibrium, excess demand for 
security raises its price, and excess supply lowers its price. In this situation, it is argued that the 
sign for rate of change of price with respect to time depends on the sign of excess demand. If the 
demand and supply functions are made linear at constant equilibrium price, deterministic model 
of price adjustment is realized with respective sensitivities. In the analysis of the solution of the 
deterministic model, it was observed that in the long run asset price settles at a constant steady 
state where no further change can occur. 

In [4] a natural decay equation is developed. The equation describes a phenomenon where a 
quantity gradually decreases to zero. In the work, it is emphasized that convergence depends the 
sign of the proportionality parameter. If the proportionality parameter is negative then it turns 
into a growth equation but if it is positive, it stabilizes in the long run. In the study of slope fields 
for autonomous equations qualitative properties of decay equation is demonstrated. It was found 
that the solution could be positive, negative or zero. In all the three cases, the solution 
approaches zero in limit as time approaches infinity. 
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In [1], dynamics of market prices are studied. It was found out that if the initial price function 
lies off the equilibrium point, in the long run stability will be realized at the equilibrium position. 
It also brings out clear case by case analysis of the solution of second order price equation by 
introducing unrestricted parameters that brings speculations in behavior of buyers and sellers. 
The author argued that depending on the signs of these parameters, buyers will cut or increase 
their purchases. Similarly, sellers will cut or increase their supplies. It was also noted that 
depending on signs of parameters, price function stabilizes in the long run. 



In [2], equilibrium solutions about a special class of static solutions are discussed. The study 
found that if a system starts exactly at equilibrium position, then it remains there forever. The 
study further found that in real systems, small disturbances often a rise which moves a system 
away from equilibrium state. Such disturbances, regardless of their origin give rise to initial 
conditions which do not coincide with equilibrium condition. If the system is not at equilibrium 
point, then some of its derivatives will be non zero and the system therefore exhibits a dynamic 
behavior, which can be monitored by watching orbits in its phase space. 



The resistance-inductance electric current circuit for constant electromotive force is modeled into 
a differential equation in [3]. The stability of the solution is studied in the long run and is found to 
be a constant, which is the ratio of constant electromotive force to resistance. It was also found 
that if electromotive force is periodic, in the long run, current executes harmonic oscillations. In 
this case, steady state solution is the periodic part of the solution. The resistance-capacitor 
electric current circuit equation is also discussed for constant electromotive force. The solution 
was an exponential function, which converges to zero in the long run. Also, in resistance- 
inductance-capacitance series circuit, a second order equation was developed. Its solution 
consists of an exponential homogeneous part and a periodic integral part. The study found out 
that the homogeneous part converges to zero as time approaches infinity, while the periodic part 
exhibits practically harmonic oscillations. 



The literature is silent and it is worthwhile developing a second order wage equation, solve it and 
discuss its stability. The solutions of linear second order ordinary differential equations are 
presented in [5; 8] using the method of variation of parameters. 
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2. Modeling second order wage equation 

In this section, we consider modeling a second order linear differential wage equation. In this 
case, the number demanded and supplied of labor is taken as functions both current wage rate 
and wage trend prevailing in the market at that time. This is because wage trend will guide 
laborers expectations regarding future wage level. The expectations will also influence labor 
demand and supply for future decision making. To include wage trend and expectations, we 
include continuous time derivatives in the model, i.e. 

™ (2.1) 
dt 

and 

^ (2.2) 
dt 

The derivative (2.1) shows whether the wage function is rising, and the derivative (2.2) shows 
whether it is rising at an increasing rate. If we consider wage trends, derivatives as additional 
arguments to labor demand and supply are taken into account. The demand and supply functions 
are thus given as 



and 



dt dt 1 

respectively. Suppose the demand and supply functions of labor are linear, then they are 

dW d 2 W / \ 
N f =n-oW + T—— + y/—^-, (n,(T>0) (2.5) 
d dt Y dt 2 V ' 

and 

N^=-0 + AW + &— + 4^-, (6U>0) (2.6) 
dt dt 

respectively. The other parameters r, y/, 3 and % introduced to dictate employers and 
laborers expectations have their signs operating in free range. For example, ifr>0, a rising 
wage rate causes the number of laborers demanded to increase. This suggest that employers 
expect rising wage rate to continue to rise and prefer to increase employment now when the 
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wage rate is still relatively low. On the other hand, for r < , the wage trend will be falling, and 
employers will opt to cut current employment while waiting for wage rate to fall even further. 
The inclusion of the parameter y/ makes employers behavior to depend also on the rate of 
dW 

change of wage rate . Thus introducing new parameters r and y/ injects wage rate 

dt 

speculation in the model. Similarly, if we let 9 > , then a rising wage rate causes the number of 
laborers supplied to fall and thus laborers expect the rising wage rate to continue rising and 
prefer to withhold their services now while waiting for the higher wage rate. On the other hand, 
for 9 < , wage rate will shows a falling trend and laborers will prefer to offer their services now 
while the wage rate is still relatively high hoping that any delay will see the wage rate fall even 
further. Therefore, introduction of the parameter £ makes the behavior of laborers to depend 
much on the wage rate just as the employers. 

We now resort to demonstrating an implicit model by assuming that only labor demand function 
contains the wage expectations. Specifically, we set both S and £, of function (2.6) equal to zero; 
while t and y/ of function (2.5) are set as non zero. Further, we assume that labor market clears 
at every point in time. We then set the number of laborers demanded equal to the number of 
laborers supplied to obtain equation 

d 2 W dW „ r , 

— - + a, + a ? W=b (2.7) 

dt dt 



with a, = — , a 2 = -( — + - j, and b = -[ - + - j , and this is the required second order linear 

¥ V ¥ ) V ¥ ) 

ordinary differential wage equation. 

3. Solution of the differential equation 

In this section, we demonstrate the solution of the differential equation (2.7) using the method of 
variation of parameters. We first solve equation (2.7) for a complementary solution by 
considering its homogeneous equation 
d 2 W dW 

a ^- + n— + a 2 W=Q. (3.1) 

dt 2 1 dt 2 

The characteristic equation of equation (3.1) is 
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r 2 +a l r + a 2 =0 . (3.2) 
Equation (3.2) is quadratic in nature and can be solved for r using the quadratic formula. The 
roots are therefore given as 

2 

If we substitute back the values of a x and a 2 as given in equation (2.7) the roots become 



(3.4) 



2^ 

The roots (3.4) yields varied results depending on the value of the 

(t)~ Ja + A'] 



discriminant 



Case I: Suppose — 



. This is discussed in three cases. 



> Othen roots r x and r 2 of the characteristic equation (3.2) are 



real and different. The complementary solution therefore becomes 
W c (t) = Cj exp r x t + c 2 exp r 2 t 



(3.5) 



T 1 + S a + ' 
WJ I ¥ 



Case II: Suppose j 
(3.2) are real and equal; that is, 

T 

r u r 7 = a = 

2 ¥ 

and the complementary solution then becomes 
W c (t) = (c l +c 2 t)expat. 



- Othen the roots r x andr 2 of the characteristic equation 



(3.6) 



(3.7) 



Case III: Suppose 



— j +4| a + ^ j < then the roots r x and r 2 of the characteristic equation 



(3.2) are complex; that is, 

r x ,r 2 =a + if3 , with a = 
and the complementary solution therefore becomes 



2 ¥ 



(3.8) 
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W c (t) = expat(c l cos fit + c 2 sin fit) (3.9) 
To solve equation (2.7) completely requires finding an integral function W,(t) by the method of 
variation of parameters. This is also discussed in three cases. 



Case 



I: Suppose f — 1 +4\ a + ^ | > the roots r x and r 2 are real and different, from 

\¥J I ¥ J 



the 



complementary solution (3.5), the integral function takes the form 

Wj (t) = Vj exp r v t + v 2 exp r 2 t 
where v x and v 2 are functions of t . The system of equations to be solved is 

v[ exp rf + v' 2 exp r 2 t = 
VjVj exp r x t + v' 2 r 2 exp r 2 t = b 

The system of equations is solved by Cramer's rule for v x and v 2 to obtain 

b 



j exp(- r x t), rj ^ 0, r x ^ r 2 



and 



V 2 = — — — | eX p(~ r \ ^ r 2' r 2 * 



Substituting solutions (3.12) and (3.13) in integral solution (3.10) gives 



W I (t) = 



r x ,r 2 *0 



From solution (3.4), we have r x r 2 = -( — + - I ; and sinceZ? 

I ¥ 

(3.14) becomes 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



the integral solution 



(3.15) 



WJt)=W =^-?-, o*-k 
a + A 

If we consider complementary solution (3.5) and integral solution (3.15), the general solution 
becomes 



W(t) = q exp r x t + c 2 exp r 2 t + W 



(3.16) 
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The general solution (3.16) is solved for a particular solution if we use the initial conditions. 

Suppose W(t)\ _„ = W a and ^—\ = then the particular solution becomes 

dt \ t=0 

W(t) = -^iw -Wlexprj —(w o -w)expr 2 t + W (3.17) 

r 2 - r x r 2~ r \ 

II: Suppose ( — 1 + 4| a + ^ \ = the roots n and r 2 of the characteristic equation are real 

W) V ¥ J 



Case 



and equal. From solution (3.7), the integral function takes the form 

Wj (t) = Vj exp(af ) + v 2 t exp(at) 
where v l and v 2 are functions of t . The system of equations to be solved is therefore 

vj exp(at ) + v' 2 t exp(ctf ) = i 

v[a exp(at ) + v' 2 (l + at )exp(at) = b\ 
This is easily solved by Cramer's rule for v x and v 2 to obtain 

Vl = (^ + ^) eXP( "^ 



(3.18) 



(3.19) 



(3.20) 



2 = exp(-af) 



a 



(3.21) 



The integral function is found by substituting solutions (3.20) and (3.21) in solution (3.18) to 
obtain 

b 



W I (t) = 

u 

But since fl.T + 1^1 = 0, f^T=« 2 =-f^ 

WJ \ ¥ ) K 2 ¥j I V 

equation (2.7), the integral function becomes 
a + A 

The general wage function is therefore found to be 
W(f) = (q + c 2 t)exp(at) + W 



(3.22) 



and because/? = 



(3.23) 



(3.24) 
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If we use initial conditions W(0| /=0 = W and 



dW \ 

~d7l 



Case III: Suppose 



W(t) = (W - Wjl - at)exp(at)+W 



= then the particular solution becomes 



(3.25) 



4| 1 |<o the roots r x and r 2 of the characteristic equation (3.2) 
W 



are complex. The roots (3.8) and the complementary solution (3.9), gives an integral function 

W, (0 = ^ exp(a + ij3)t + v 2 exp(a - ij3)t (3.26) 

where v t and v 2 are functions of t . We use Cramer's rule to solve for v t a«<i v 2 in the system 
v[ Qxp(a + ifl)t + v' 2 exp(a-i/3)t = (a)] 



^'(a + i/?)exp(a + + v' 2 (a - z'/3)exp(a - ifi)t = b (b)\ 



3 - ia) 



exp(- (a + //?>) 



^exp(- (a -i/3)t) 



Substituting solutions (3.28) and (3.29) in the integral function (3.26) gives 
b 



W,(r) = 



a 2 +j3 2 



(3.27) 



(3.28) 



(3.29) 



(3.30) 



But from solution (3.8), a 2 + (5 2 = -\ a + ^ j and sinceZ? = -( ^±IL j the integral function (3.30) 
\ ¥ ) V ¥ ) 



Wj(t)=W 



e+T] 



a * -A. 



(3.31) 



a + A 

The most general solution is thus written as 

W(t) = expctffo cos fit + c 2 sin J3t)+W. (3.32) 
Solution (3.32) is solved for a particular solution if we use the initial conditions. Suppose we let 
dW\ 

W(t)\ t=Q = W and — — = the particular solution becomes 
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W(0 = exp(ctfl (W -W) cos fit (W Q -W)sin/3t \ + W 



(3.33) 



4. Results, analysis and interpretation 

In this paper, a second order wage equation 
d 2 W dW 



dt 2 



dt 



- + a ? W = b 



(4.1) 



with a, - 



-m 



andb = 



has been developed for the first time. In this 



¥ V ¥ ) \ ¥ J 

equation, > is a parameter that shows the number of laborers supplied that does not depend 
on the wage rate, X > is a parameter that shows the proportion by which the number of laborers 
supplied responds to variation in wage rate, 77 > a parameter that shows the number of laborers 
demanded that does not depend on wage rate, and <r>0 a parameter that shows the proportion 
by which the number of laborers demanded responds to variation in wage rate. The parameters 
r and \j/ introduced to dictate employers' expectations have their signs operating in a free range. 
For example, if r > , a rising wage rate causes number of laborers demanded to increase. This 
suggests that employer expects a rising wage rate to continue to rise and prefer to increase 
employment now when the wage rate is still relatively low. On the other hand forr <0 , the wage 
rate shows a falling trend and employer will opt to cut employment now while waiting wage rate 
to fall even further. The parameter y/ makes the employer's behavior to depend on the rate of 
change of wage rate. 



Equation (4.1) has been solved in this paper by the method of variation of parameters for the first 

time. When | — ] + 4\ cr + ^ ] > , the general solution was found as 

¥) \ ¥ J 



W(t) = c l expr l t + c 2 expr 2 t + W (4.2) 
with r x and r 2 described in solution (3.4), and W in solution (3.15). In this case, 



suppose i// > 0, then A\ a + ^ j > and ( — j +4\ a + ^ 1 > 0, V values of 
V ¥ J WJ I ¥ J 



t . The solution (4.2) is 
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therefore valid. Moreover, with y/ > 0, since <r, A > , ( — 1 + 4\ <J + ^ j i s positive and its square 

WJ K ¥ J 

root exceeds — . Therefore solutions (3.4) produce one positive root r. and one negative 

WJ 

root r 2 . Consequently, inter temporal equilibrium is dynamically unstable. For the function (4.2) 
to be stable, we have to set constant q to zero so that it becomes 

W(t) = c 2 expr 2 t+W (4.3) 
dW\ 

If we use the initial conditions W(0| t=0 = W and — — | = then solution (4.3) becomes 

W(t) = ^V -w)expr 2 t + W (4.4) 
The solution (4.4) is now investigated for stability by taking the limit as t tends to infinity, i.e. 

W(t) = lim, Jw ~w)expr 2 t + w\ (4.5) 
In this case, (w o - w)is constant and the value of limit function (4.5) depends on the exponential 
factor exp r 2 t . In view of the fact that r 2 < 0, (w o - VV^exp r 2 t — » as t —> oo . The limit function 
(4.5) therefore becomes 

W(t)=W. (4.6) 
This means time path of the wage function (4.4) consequently moves towards equilibrium 
position in the long run. 

Considering particular solution (3.17), suppose we let^ < with a, A >0, the expression under 

the square root in solution (3.4) is less than — and the square root must be less than — ; 

WJ ¥ 

therefore for ( — 1 + 4\ G + - 1 > 0, if we let z < 0, then the roots of (3.4) would produce two 

WJ K ¥ J 

negative roots. The solution (3.17) is therefore investigated for stability by finding its path, i.e. 
W(t) = lim,^ j^— ^— j(w - w)exp r x t - lim,^ j^— ^— j(w o - w)exp r 2 t + lim,^ W. (4.7) 
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In this case, — - — IW -W)and — - — IW -VK)are constants and the value of the limit 

function (4.7) depends on the exponential factors exprj and expr 2 t. In view of the fact that 
r p r 2 < Othe limits of the first and second term of function (4.7) both tends to zero; thus 

W(t)=W. (4.8) 
This means that the wage function (3.17) consequently moves towards the equilibrium position 
in the long run and it is therefore dynamically stable so long as y/ < and r < . 

Interestingly, if we consider f — 1 + 4| — + - 1 = , the general solution of equation (4.1) is as 

WJ I v J 

shown in (3.7). If we use the initial conditions, a particular solution (3.25) is obtained and it is 
investigated for stability, by taking its limit as time tends to infinity i.e. 

W(0 = lim, Jw -w\l-at)exp(at)+w\ (4.9) 
In this case, the first term of the limit function (4.9) consists of a linear 
factor(w o -W^L- at) and an exponential factor exp(at). Its value therefore depends on the 
exponential factor exp(ctf ) . In view of the fact that a < 0, (w o -W^l- at) exp(at) — » as t — > oo . 
The limit function (4.9) therefore becomes 

W(t)=W. (4.10) 
This means the time path of the wage function (3.25) consequently moves towards the 
equilibrium wage rate as time tends to infinity. 

Further analyses of the limit functions (4.5), (4.7) and (4.9) is possible by considering the 
relative positions W Q andW ; that is, by comparing the relative positions of the initial wage rate 
and the equilibrium wage rate. This is discussed in three different cases. 

CASEI: In this case, we consider both the limiting functions and let W = W .This means the 
limit functions (4.5), (4.7) and (4.9) becomes W(t) = W at infinite time, which is a constant path 
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and is parallel to the time axis. The wage function in both situations becomes stable at 
equilibrium wage rate in the long run. 

CASE II: In this case, we let W > W . The first term on the right hand side of function (4.5) is 
positive but it decreases since as t — > qo it is lowered by the value of the exponential factor 
expr 2 ? forr 2 < 0. The first term on the right hand side of function (4.7) is positive if r, > r 2 and 
the second term is only positive if r x < r 2 .Therefore, they will decrease since as t — > co they are 
lowered by the values of the exponential factors expr^and exp /^respectively. Finally, the first 
term on the right hand side of function (4.9) is positive but it decreases as t — > oo since it is 
lowered by the exponential factor exp(ctf). The limit functions (4.5), (4.7) and (4.9) thus have 
their time path asymptotically approaching the equilibrium wage rate W from above, and in the 
long run, become stable. 

CASE III: In this case, we let, W < W i.e. the initial wage rate is taken to be less than the 
equilibrium wage rate. The first term on the right hand side of the limit function (4.5) is negative 
and the exponential factor infinitely makes W Q to rise asymptotically towards the equilibrium 

wageW ast — >ao . Similarly, the first term on the right hand side of function (4.7) is negative if 
r x > r 2 and the second term is only negative if r x < r 2 , making W to rise asymptotically towards 
the equilibrium wage rate W as t — > qo. Finally, the first term on the right hand side of function 
(4.9) is negative and it infinitely makes W to rise asymptotically towards the equilibrium wage 
W as t — > oo . These three cases are illustrated in figure 4.1. 
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W(r) 



^0 




Figure 4.1: Stability Analysis of the Wage Function 



Figure 4.1 shows that when W = W, then W(t) = W , which is a constant function. 
If W > W, thenW(t) decreases asymptotically towards W, and if W < W thenW(t) increases 
asymptotically towards the equilibrium wage rate W. The results therefore shows that as t — > <x>, 
the functions (3.18) and (3.26) approach the equilibrium wage rate and becomes stable so long as 
y/ < and r < . 

We now turn to investigating the solution when [ — 1 + 4| °" + - ] < . In this case, the general 
wage function (3.32) has been developed. If the initial conditions are used, a particular function 

(3.33) is obtained. The function (3.33) is now investigated for stability, i.e. since a = — — is the 

2y/ 

real part of the complex root, if we \ety/<0 and r <0, thena<0. The wage function (3.33) 
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therefore becomes dynamically stable. The time path in this case is one with periodic fluctuation 
of period ^ ; that is, there is a complete cycle every time t increases by ^ , where /? is as 

defined in (3.8). In view of the multiplicative factor expctf the fluctuation is damped. The time 
path, which starts from the initial wage, W(f)| t=0 = W converges to an inter-temporal 
equilibrium wage W(t) = W in a cyclical fashion. This is illustrated in figure 4.2. 

way 




W(t) = W 



Figure 4.2: Periodic Stability Analysis of the Wage Function 



Figure 4.2 show that sinceor<0, the exponential factor expo* continually decreases as 
t — > go and each successive cycle gives smaller amplitude than the preceding one and the ripples 
naturally dies slowly. 

4. Conclusion 

In this paper a second order wage equation has been developed and solved using the method of 
variation of parameters for the first time. The subsequent wage functions have been analyzed and 
interpreted for stability. The equation incorporates speculative parameters operating in a free 
range in its modeling. The parameters are introduced to dictate employers' expectations on future 
wage rates. The variations of these parameters have caused stability and instability in the wage 
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function in certain circumstances. Above all, where the wage function takes an exponential form 
with particular assumptions, as time approaches infinity it asymptotically approaches the 
equilibrium wage rate. Where as in a case of an exponential and a periodic factor, the time path 
shows a periodic function whose successive cycles elicit smaller amplitudes and the ripples 
eventually dies naturally at equilibrium wage rate as time approaches infinity. It has been 
realized that where the wage rate is determined by free market forces of demand and supply, 
volatility in wage rate may be observed if not controlled. This may increase uncertainties and 
cause anxiety about investment and employment in the economy. The paper therefore proposes 
government intervention by creating a middle path in which wage rate is allowed to oscillate 
freely within a narrow band managed by employers in consultation with the workers under the 
watch of the government. 
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